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Abstract 

Matrix pair beamformer (MPB) is a blind beamformer. It exploits the temporal structure of the signal 
of interest (SOI) and applies generalized eigen-decomposition to a covariance matrix pair. Unlike other 
blind algorithms, it only uses the second order statistics. A key assumption in the previous work is that 
the two matrices have the same interference statistics. However, this assumption may be invalid in the 
presence of multipath propagations or certain "smart" jammers, and we call it as matrix mismatch. This 
paper analyzes the performance of MPB with matrix mismatch. First, we propose a general framework 
that covers the existing schemes. Then, we derive its normalized output SINR. It reveals that the matrix 
mismatch leads to a threshold effect caused by "steering vector competition". Second, using matrix 
perturbation theory, we find that, if there are generalized eigenvalues that are infinite, the threshold 
will increase unboundedly with the interference power. This is highly probable when there are multiple 
periodical interferers. Finally, we present simulation results to verify our analysis. 
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I. Introduction 

Beamforming is a spatial filter, which combines the outputs of multiple sensor elements by an appro- 
priately designed weight vector so as to pass a signal of interest (SOI) while rejecting interfering signals. 
Since the pioneering work of Howells [1], Applebaum [2] and Widrow [3], it has been intensively studied 
during the past decades and widely applied in radar, sonar and wireless communications [3]-[ll] etc. 
For a comprehensive review, we refer to [12]— [15] and the references therein. 

There are various forms of implementations for a beamformer. Some use the direction of arrival (DOA) 
of the SOI and directly calculate the weight vector by sample matrix inversion (SMI) [16]. Some employ 
a reference signal (e.g. training signal [17] and decision feedback signal [6]) to iteratively calculate the 
weight vector. And there are also blind beamformers which do not require the DOA or the reference 
signal [18]-[31]. 

To achieve blind beamforming, we need to exploit the special properties of the SOI, such as constant 
modulus, non-Gaussianness, high order statistics etc. One important property is the temporal structure 
of the SOI. Such kind of blind beamformers are extensively studied for CDMA systems, where the 
inherent structure of the spreading codes can be exploited [23]-[26], [28]— [31]. The advantage of this 
scheme is that it only relies on the second order statistics of the covariance matrix pair. Although the 
implementation details differ, the main idea of these approaches is to exploit a pair of array covariance 
matrix and hence will be referred to as matrix pair beamformer (MPB) in this paper. The MPB projects 
the discrete sequence in each antenna onto the space spanned by the SOI's signature vector (for a CDMA 
system, it is the spreading code vector) and another carefully designed base vector. Then two sets of array 
snapshots (i.e. signal snapshot and interference snapshot) are acquired to calculate a pair of covariance 
matrices. With the processing gain, the desired signal power in the signal snapshot is generally greater 
than that in the interference snapshot. A key and common assumption is that the interference statistics 
in the two snapshots are identical. These two features enable the separation of the signal steering vector 
and the interference covariance matrix from the two matices. And the weight vector derived from the 
dominant eigenvector of the matrix pair will maximize the output signal to interference plus noise ratio 
(SINR). 

However, this key assumption that the two matrices share the same interference statistics is not valid in 
many cases, which we refer to as matrix mismatch. Matrix mismatch may occur when there are interferers 
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with certain periodical structure, such as multiple access interference (MAI) in CDMA systems, tones 
and some other "smart jammers" in radar systems and it may lead to the failure of a system. In practical 
applications, to avoid such an unexpected failure, it is necessary to understand the detailed performance 
of a scheme when the assumption/condition is not satisfied. To our best knowledge, little effort has been 
devoted to analyzing the effect of matrix mismatch on the performance of MPB. This paper aims to 
analyze the MPB's performance under matrix mismatch, and our contributions are: 

• proposing a general framework to model various existing MPB schemes; 

• deriving analytical expressions for the normalized output SINR as the performance measure; 

• discovering a threshold effect for MPB, i.e., due to matrix mismatch, the performance of MPB 
degrades rapidly when the input signal to noise ratio (SNR) is below a predicted threshold, and the 
main beam points to the directions of the interferers; 

• explaining how MPB works "blindly" by "steering vector competition"; 

• discussing various factors that have impact on the threshold, and showing that when there is an 
generalized eigenvalue that is infinite which is called the noise-free covariance matrix pair, the 
threshold SNR increases unboundedly with the interference power; 

• discussing several typical scenarios and showing that MPB is very vulnerable to multiple periodical 



The rest of the paper is organized as follows. Sec. II presents a general framework of MPB to cover 
and reinterpret the basic ideas in [23]— [31], followed by a formulation of the matrix mismatch problem. 
In Sec. HI, we present the expressions for MPB's weight vector and its normalized output SINR, which 
uncovers the inherent threshold effect caused by matrix mismatch. The discussion relies heavily on the 
approximation of the generalized eigenvalue, which is derived in Appendix A using Gerschgorin theorem 
[32], [33]. Sec. IV applies the Weyl-Lidskii type theorem in matrix perturbation theory [32] to analyze 
MPB's threshold SNR, and discusses it in two typical scenarios. Finally, Sec. V presents simulation 
results to verify our theoretical analysis and Sec. VI concludes the whole paper. 



A. Signal Model 

Consider an antenna array of L isotropic elements that receives D + 1 signals from a far field. After 
preprocessing (e.g. mixing, filtering, etc.) and sampling, the output Lxl array vector can be written as 



interferers. 



II. Problem Formulation 



D 




(1) 



i=0 
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Fig. 1. Block diagram of the matrix pair beamformer in a CDMA system. Data corresponding to the k symbol are blocked and 
projected onto two spaces to produce signal channel snapshot xs(fc) and interference channel snapshot Xi(fc), respectively. 



where Sj(n) is the discrete sequence of the ith signal with normalized power; Pi is its power; is 
the steering vector, which depends on the DOA and the array geometry; v(n) is the additive white 
Gaussian noise (AWGN) vector with zero mean and covariance matrix <r 2 I. Assume so(n) is the SOI, 
and si(n), S2 («),..., szj(n) model all the possible interferers like MAI and jammer etc. This paper 

considers the case when so(n) has some inherent temporal structure expressed as 

+00 

so(n) = ^2 b (k)c [n - r(k)}, (2) 

fc=— oo 

where co(n) is known and supported on < n < N — 1. Thus, so(n) can be viewed as a train of pulses, 
with each of them being a delayed and scaled replica of co(ra). And r(k) and bo(k) are the corresponding 
delay and amplitude of the kth pulse. This model is common in modern communication systems. For 
instance, in CDMA system, co(n) is the spreading code of the desired user, bo(k) is its data bits and r(k) 
is the delay of the kth symbol. Without loss of generality, let us consider bo(k) 6 {±1}, co(n) 6 {±1} 
and r(k) = kN + uq. Then N is the processing gain and no is the propagation delay. 

B. A General Framework of the Matrix Pair Beamformer 

We first propose a framework called matrix pair beamformer (MPB) to cover the common ideas in 
[23]— [3 1]. Our strategy is to use orthogonal projection operation to model their ways of estimating the 
covariance matrix pair. By this framework, we will find a common threshold effect in these methods. 
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The steering vector a^ in (1) is a spatial signature of the zth signal, which is different from other a.j so 
long as they arrive from different directions. Beamformer is a spatial filter that exploits such difference 
to pass the desired signal So(n) while suppressing s\(n) . . . S£>(n) and v(n). A statistically optimum 
beamformer [12], [15] generally requires at least, either explicitly or implicitly, the information about 
the steering vector ao and the interference covariance matrix. The latter one may be replaced by the 
data covariance matrix, so the remaining problem is how to acquire ao- In DOA-based beamformer, it is 
calculated by the DOA and array manifold information. As for training-based method or decision directed 
method, it is inherent in the correlation vector between the reference signal and the data vector. 

To work "blindly", i.e. without any explicit information of DOA, the methods in [23]— [3 1] exploit 
the SOI's temporal signature co(n) to acquire these spatial statistical information. Specifically, it is 
implemented by two orthogonal projections and a generalized eigen-decomposition of a covariance matrix 
pair. Hence, we refer to them as matrix pair beamformer (MPB) in this literature. In Fig. 1, we summarize 
the common structure of MPB. With the data segmentation, the array outputs corresponding to the kth 
symbol of the desired user can be expressed in the following matrix form: 

/Pbbo(kj\ aoc^+A/efsf (k)+V(k) 



X(fc): 



x(kN + n ) 



x(kN + n + N - 1) 



(3) 



Z(fc) 



where 



A/ 

co 
Si(k) 



ai a 2 • • • a D 

co(0) c (l) •• 
si(fe) s 2 (k) ■ ■ 



Si(k)= Si(kN + n ) 
V(k) = \ v(kN + n ) 



co(N-l) 
s D (k) 

• Si(kN + n + N - 1) 

• v(kN + n + N - 1) 



ej4diag{Pi,iV-- ,Pd}. 

Then, the fcth data block in each antenna is projected onto two subspaces: signal space S spanned by the 
SOI's temporal signature vector hs = co/VN, and a specifically designed interference space Z=1Z{Kz}, 
respectively. Without loss of generality, assume the columns of HiGC Arxrx are orthonormal. Then, the 
projections produce the signal snapshot x.s(k)=X.(k)hg and the interference snapshot Xj(fc)=X(A;)HJ. 
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Define <jg o =NPo and ct| o =^-||Hj co|| 2 . Then, the covariance matrices of x$(k) and Xj(fc) are 

Rs 4 E{xs(A;)xf (k)} = 4o a oa? + Qs (4) 

Rx 4 1e{X x (A:)X^ (fc)} = <7£aoa£ + Qx (5) 

where cr| and af are the SOI's powers in the signal channel and interference channel (c.f. Fig. 1), 
respectively. Qs and Qx are the covariance matrices of the interference-plus-noise in them, defined by 

Qs = E{Z(k)h* s h^Z H (k)} = A,* 5 Af + a 2 I, (6) 

Qx = E{Z(k)U* s li^Z H (k)} = Aj^xAf + a\ (7) 

where *s = cx 2 INR • $ So , * x = cj 2 INR • * Io , INR = Pi/cr 2 , and 

<l>s = n)E {S?(fc)h$h£SK*0} n), *x„ = ^«|E {S?(fc)H£H£SKA:)} «), (8) 

with fi/ = (e//o- 2 )/INR. We can see that flj is independent of INR and only depends on the relative 
strength of the interferers. So are <l?s and <&x > anci we will use this conclusion in Sec. IV. 

The MPB uses the eigenvector corresponding to the largest generalized eigenvalue of the matrix pair 
(Rs, Rr) as the weight vector w, which is the solution to the following equation 

R 5 w = A max R x w, (9) 

where A max is the largest generalized eigenvalue of (Rs, Rr). w is applied to xs(A;) to yield the output 

y (k)= W H ^ s (k)=[a s Mk)} ^ H a^ H Z(k)h% (10) 

" v ' " v ' 

Vs(k) yi(k) 

where ys(k) and yi(k) are the output signal and interference-plus-noise. Eq. (9) is also equivalent to 

(Rs - Rr) w = (A max - l)R x w. (11) 

Using (4) and (5), we have 

Rs - Rx = (<4 - 4 o )a a^ + Qs - Q x . (12) 
Since the columns of Hj are orthonormal, i.e. Hj Hj = I, its spectral norm is one and we can have 

2 A -Pomxt// ,|2 ^ ^0||tt ,|2|, |,2 ^0 ^ 2 

rx r x r x 

The above inequality holds strictly when choosing Hj / hs (either rx > 1 or Hj / hs for rx = 1). 
Furthermore, it is commonly assumed [23]— [31] that Qs = Qx- Then, Qs~ Qx = O in (12), and the 
dominant eigenvector of (Rs — Rx, Rx) = ((f"s ~~ °i ) a o a ^^ ^) * s 

w opt = mR^^o = //Q^ao, (13) 
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where \i and // are scalars. Then w opt will maximize the output interference plus noise ratio (SINR) 
[15], and the optimal SINR is 

£{\ys(k)\ 2 } 



SINR C 



al^Q^ao (14) 



opt £{\yi(k)\ 2 } 

All the methods in [23]— [3 1] share the structures described above. They only differ in the dominant 
eigenvector searching algorithm and the interference space X (i.e. Hi). In most existing approaches, 
X is a one dimension space (rj = 1). The pre- and post-correlation (PAPC) scheme [23]-[27] uses 
Ri = E{x(n)x ff (w)}, thus it is equivalent to choosing Hj as 



H 7 



••• 1 ••• 



T 



(15) 



where only one component in Hj is nonzero. The Maximin scheme in [28] and [29] employs a monitor 
filter to isolate the interference, which can be interpreted as 



H x = c 



I e j27r/ M F . . . e j2irf M F(N~l) 



T 



(16) 



where /mf £ (0, 1] is the normalized center frequency of monitor filter, and is the Hadamard product. 
C. Matrix Mismatch 

We see that the MPB relies heavily on the key assumption that Q$ = Qj, namely the interferers have 
the same second order statistics in the two channels. By (6)-(8), we know it is valid when each interferer 
is random enough in the temporal domain, say directional white noise (c.f. Sec. IV-B1). However, it is 
generally not satisfied, especially when there are multiple deterministic periodical interferers, like tones 
and MAI (c.f. Sec. IV-B2). This is because when the interferers are deterministic and periodical, the 
expectations in (8) can be eliminated. Then Q$ = Qj requires S^{k)h* s = Sj (k)H x , which is highly 
improbable when hs^Hz- We term this as "matrix mismatch". To our best knowledge, very little effort 
has been devoted to analyze this problem. Therefore, we will investigate the performance of MPB in this 
more general case. Before we proceed, we define the normalized output SINR to measure performance 
degradation with respect to that of no matrix mismatch. 

Definition 1: The normalized output SINR is defined as the actual output SINR of MPB normalized 
by the optimal value, i.e. 

1 



q A £{\ys(k)\ 2 } 



E{|y/(fc)| 2 } 



SINR 



opt 



(17) 



where SINR opt is given by (14), and w is the solution to (11) without the assumption of Q$ = Qj. 
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G generally depends on the input SNR and the interference powers. So a reasonable way to characterize 
the performance is to plot G against the input SNR = a^ o /a 2 , when fixing the interference powers. In 
the following sections, we will base our analysis on G(SNR), which we will refer to as operating curve. 
Moreover, we assume infinite sample size so that the finite sample effect is ignored. 

III. Performance Analysis of the Matrix Pair Beamformer with Matrix Mismatch 
In this section, we will derive the operating curve of MPB, and discuss how it works blindly. 



A. Operating Curve of Matrix Pair Beamformer 

We base our discussions on the following assumptions, and summarize the main result in theorem 1. 

Assumption 1: The spacing of DOA between any two signals is large enough (greater than a mainlobe), 
so that {ai}fL are linearly independent and the projection of ao onto span{aj}^ =1 is much less than 
ll a o||- 

Assumption 2: The steering vectors of all signals are normalized so that ||aj|| 2 = L, (i = 0, 1, . . . , D). 
Theorem 1 (operating curve): The normalized output SINR of MPB with matrix mismatch is 

G(SNR) (tw^snrFTI - SNR>SNRt2 (18) 

. [l-SNR/SNR T o+^o°(L/3SNR/iV+l)J Gl ' SNR < SNR ti 
where G\j and Gl are the normalized output SINR when SNR=+oo and SNR=0 (— oodB), respectively. 
P—Na^/ag is the normalized power leakage ratio (PLR) in interference channel. Pt=1/Gjj—\ is the 
output interference to noise ratio, and SNRjo is the empirical threshold SNR. Their expressions are 

a^Q/ao a^Q/QsQ/ao 
/3+$/SNR TO ^ N-fj 



K °= N-0 v 1 ~fT~ J <21) 

where (•)+ = max{-,0}, and 71 is the largest nonzero generalized eigenvalue of (Q5— Qx, Qz). 1 The 

SNR at which G(SNR) is close to Gl and Gjj (within 3dB), which are given by 

-1 




SNRti =(1-i/^| SNRto, SNR T2 = I 1 - \l I SNR T0 , (22) 



2P / + 1 / 

'if (Qs Qi, Qi) has less than D nonzero eigenvalues, pad them with zeros up to D and order them decreasingly. 
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respectively. 

Though we will give expression for Gl in Sec. III-D, its specific values are of no interest to us. Fig. 
2(a) shows a typical curve of G(SNR). We plot the curve in failure area and operating area given in 
(18), and connect their ends by a dashed line. We can see that the performance of beamformer degrades 
rapidly when the input SNR is below SNRj2- And it fails completely after reaching SNRji- Therefore, 
matrix mismatch causes a threshold effect in MPB, and SNRj2 is a critical parameter to be optimized. 

There are two special cases of G(SNR). The first one is SNRti=SNRt2=0 (— oodB). This happens 
when there is no matrix mismatch, i.e. Q,s=Qi, so that 71=0 and SNRjj given by (20) and (22) are zero. 
Then, (18) implies the operating curve is a horizontal line as shown in Fig. 2(b). Another interesting case 
happens in PAPC schemes mentioned in II-B, whose PLR is /3=1>0. If it further satisfies (N — j3)/ji<(3, 
then by (20) and (22), SNRt;l=SNRt2 = +°o and G(SNR) only has the failure area, which decreases 
in the order of 0(SNR~ 2 ) (c.f. Fig. 2(c)). In the rest of the section, we will derive G(SNR) and reveal 
how MPB works blindly under matrix mismatch. The discussion of SNRxo is left to Sec. IV. 

B. Derivation of the Weight Vector for MPB 

As a step to prove Theorem 1, we first derive the expression of the weight vector for MPB. By the 
arguments in Sec. II-B, it is the solution to (11). Define *a — *s — *x- Then, by (6), (7) and (12), the 
problem becomes solving the following generalized eigenequation 

[(4 - 4 )*o*o + A/* A Af ] w = (A max - l)R x w, (23) 

where A max is the largest eigenvalue of (Rs, Rj). Sec. II-B already gave the result for w when there is no 
matrix mismatch, i.e. 3>a=0 (or Q,s=Qx). We concluded that it is optimal in the sense of maximizing 
the output SINR. However, in the presence of matrix mismatch, we have <J>a/0, which is the key 
challenge for solving (23). To deal with this problem, we first have the following two observations: 
• By left-multiplying A 1 _ 1 Rj 1 to both sides of (23), we can see that w can be expressed as 

D 

w = ffcR^ao + ^ r/iR^a; (24) 

i=i 

i.e. it is a linear combination of R^ao, Rj 1 ai, . . . , R^a^i, where i]o = ^° a^w, a« is the ith 
column of Aj, and rji is the ith component of rj = $aA^w. To determine r]o, rji, . . . , rjo, we only 
need to substitute (24) back into (23) and solve a linear equation. However, we will immediately 
discover that the solution is intractable for further analysis because A^Rj 1 A/ and 3>a are not 
diagonal. To overcome this, we need the next observation. 
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Fig. 2. (a) A typical curve of G(SNR) has failure area, threshold area and operating area, separated by SNRti and SNRj2- 
(b) Without matrix mismatch, the curve has only operating area, (c) A curve has only failure area for some cases in PAPC. 



Suppose we can factorize 3>a into diagonal form such that A/3>aA^ = A e rA e and A^R 2 - 1 A e = 
I, where T is a diagonal matrix. Then, by repeating the above procedure except for replacing 



ax, . . . ,&d by a £l , . . . , a. tD , with a ei being the ith column of A e , we can solve w as 

D 



W 



i=l 



1-7, 



R^a, 



(25) 



where 7, is the ith diagonal component of T, r/o is an arbitrary constant, and ^ = a^R^an. 
So far, the only thing left is to figure out a way to factorize <&a so that it meets the above require- 
ment. Consider the generalized eigen-decomposition of the matrix pair (<&a> (Af^R^ 1 A/) -1 ). By the 
simultaneous diagonalization theorem (c.f. [34], pp.133), there exists a nonsingular matrix T such that 

T H * A T = r, T H (AfR^ 1 A/) _1 T = I, (26) 
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where T = diag(7i, 72, . . . , 7d) is a diagonal matrix whose diagonal terms are the generalized eigenvalues 
of (* A , (Af R^A/) -1 ). Define A e = A^T -1 )^ Then, by (26), we can verify that 

A 7 * A Af = A / (T~ 1 ) H rT" 1 Af = A e rAf 

A £ R^ 1 A e = T _1 Af Rf ^/(T- 1 )^ = I. 



Therefore, (26) is the exact decomposition of <1?a we are looking for. Furthermore, by A e = A/(T 



-1\H 



we notice that each column of A e is in fact a linear combination of ai, . . . , &£>. This, together with (25), 
implies that w is still a a linear combination of R^ao, Rj 1 ai, . . . , R^a^, just as that in (24). 



C. How MPB Works Blindly 

To fully understand the behavior of MPB given by (25) and how it works blindly, we still need the 
expressions of A max and 7$. However, A max and 7$ are the solutions to 

det {(A - 1)R X - (a 2 So - crf o )a a^ - A 7 * A Af } = (27) 
det{A(AfR^ 1 A / )- 1 - * A } = (28) 

which are polynomial equations of degree D + 1 and D, respectively. It is known that there are no 
general closed-form solutions if their degrees are higher than four. Thus, our approach here is to derive 
approximate expressions for A max and 7$. The main idea is to transform (27) and (28) into eigenvalue 
problems of a diagonal matrix perturbed by a small term. Then we can approximate the eigenvalues by 
the diagonal entries, and bound the error using matrix perturbation theory. 

We first discuss A max . Before this, let's introduce the following identity which will be used repeatedly 
in the analysis and can be derived from Properties 16 and 17 in [34, pp.5]. 



det(AI - XY) = \ m ~ n ■ det(AI - YX) 



(29) 



where X G C mxn and Y G C nxm . Substituting the factorization of $ A in (26) into (27) and using (29) 
as well as A e = A/(T _1 ) H , we can have the equivalent form of (27) as below 



(A-l 



,L-D-1 



det Rr • det < 



(A - 1)1 



4, - ° 2 







J 



AfRT 1 ' 



.A. £ 1^-7- .A. f 



0. (30) 



Then, (30) implies that the solution of (27) is the same as the eigenvalues of the following matrix 



M = 



70 + 1 7oV i>T 
Ttp T T + 1 



±£SNR+1 



(31) 
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except for the multiplicity of ones, where we have used the second equality of (26). And 70 and i^t are 
70 = (4 ~ 4>o^x^o, = (^SNR + l) Af R^a (32) 

Before proceeding on, we cite the following two lemmas. The first one is summarized from Lemma 1 
and Lemma 2 in [35], and the second one can be found in [33, pp.344]. 
Lemma 1: The quantities a^R^ao and ipT have the following results 

h -1 L N „ , F2 FT 

a ° ^ ao * ^L^SNRTiV' IM * V ^° « V ^ 

where (3 is the PLR defined in Theorem 1, and kq is a very small positive number independent of SNR. 

Lemma 2 (Gerschgorin): Let matrix A G £ nxn and a^- be its i,j th element. Define 

n 

Ri(A) = ^2 \ a i,j\i l<i<n 
as the deleted absolute row sums of A . Then all the eigenvalues of A are in the union of n disks, i.e. 

n 

X(A)e\jG l (A), 

i=i 

where Gi(A)={z<EC:\z-au\<Ri(A)} is the ith Gerschgorin disk. Moreover, if a union of k disks forms a 
connected region disjoint from all the remaining disks, then there are k eigenvalues of A in this region. 

By Lemma 1, the off-diagonal terms of each row in (31) are much smaller than the corresponding 
diagonal terms. Then, according to lemma 2, its eigenvalues are approximately 70+I, • • • ,7d+L We will 
give a bound for the approximation error in the next subsection. Now we directly use this approximation 
to discuss the behavior of MPB. Without loss of generality, assume 71 > . . . > 7d> then 

Amax ~ max{7o + 1, 71 + 1, . . . , Jd + 1, 1} = max{7o + 1, 71 + 1}. (33) 

Substituting the expression of a ¥ Rj 1 ao in Lemma 1 into its definition in (32), we can have 

T SNR 

-"• = (N -V V SUR + N <M) 
which is a monotonically increasing function of SNR. Furthermore, we will show later that 71+I is 
approximate to one of the generalized eigenvalue of (Q,s,Qx) and is almost independent of SNR. 
Hence, there is a threshold SNRjo such that A max switches from 71 + 1 to 70 + 1 when SNR exceeds 
SNRjo- Its expression is in (20) and can be derived by setting (34) to 71 and solving SNR. As a result, 
• When SNR > SNRjo, A max « 70 + 1. If 70 + 1 can be much greater than 71 + 1 as SNR increases, 
then (25) implies w ps / uiRj 1 ao- Hence, ao is the dominant steering vector, and the beamformer 
can operate properly by steering the mainlobe to the direction of the desired signal. 
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• When SNR < SNRjq, A max ~ 71 + 1. Thus w « / LtiR^ 1 a ei , and a ei is the dominant steering 
vector in (25). By A e = Aj(T~ r ) H , a ei is a linear combination of all the interference steering 
vectors. Therefore, the beamformer fails for its mainlobe pointing to the directions of the interferers. 
If then there would be SOI in Rj, which makes MPB treat it as interference and null it out. 
This explains the reason why G(SNR) decreases in the order of 0(SNR~ 2 ) in Fig. 2(c). 

Now, we further analyze how this threshold effect happens. By the factorizations (26), we can rewrite 
R5 — Rx on the left hand side of (23) as 



D 



R 5 - Rx = 70 • a So af o + ^ 7i • a e t a^f 



(35) 



i=i 



where as = ao/[a ¥ Rj 1 ao] 2 . It is obvious that 70 is the generalized eigenvalue of the matrix pair 
((c| — 0x o ) a o a cM R-z), with R^ 1 a5 being its "normalized" eigenvector. By applying (29) to (28), we can 
also verify that 71,72, ■■• ,7r> are the generalized eigenvalues of (Aj$AAf,Ri), with the eigenvectors 
being {Rj 1 a £i }^ =1 . Comparing (35) with the left hand side of (23), we can find that 70 actually measures 
the mismatch of the desired signal between the signal channel and the interference channel in Fig. 
1, and {^ji}f =1 measures that of the interferers. From the previous discussion, we know that A max = 
max{7o + 1, 71 + 1}, and w takes R^ao or R^ 1 a £l depending on which one of 70 + 1 and 71 + 1 is 
larger. Therefore, MPB blindly chooses the one with the largest mismatch metrix in (35) as its steering 
vector. We term this as "steering vector competition" . The threshold happens when the SNR is large 
enough so that 70 exceeds 71. If there is no mismatch, i.e. 3>a = O so that 73 = 0, then the desired 
signal always wins out in the competition (70 > 71), and MPB can work properly for all SNR. 



D. Proof of Theorem 1 

Let us return to our original problem of operating curve for MPB and prove Theorem 1. We first use 
the expression of w in (25) to derive G(SNR) in term of A max . The result is given by the following 
lemma, with its proof in Appendix B. 

Lemma 3: The normalized output SINR in (17) can be expressed as the following function of A max 

N 



G(SNR) 



1 + L/?SNR + AT^ (Amax) 



L/3SNR /2M 
,^ (Amax) " L/3SNR + iV^ 



+ 



1 + 



N 



L/3SNR + N 



i>S (Amax) 



(36) 



where V'S'(Amax) and V>/(A m ax) are in the following forms with ipr i being the ith component of i^t- 



a 2 ° 

Y>s(Amax)=y- Yj 



Amax 1 
L ^A m ax-(7i+ 1 ) 



, a 
'~~L 



v 



i=l 



Amax 1 
Amax-(7i+l) 



(37) 



June 7, 2011 



DRAFT 



14 



To further simplify (36), we need the expression of A max . In the previous subsection, we used the 
approximation (33) to analyze MPB. The following lemma quantifies how precise it is by bounding the 
approximation error. Besides checking the validity of (33), this is also critical in deriving G(SNR). 

Lemma 4: The largest generalized eigenvalue of (Rs,Rx) defined in (4) and (5) satisfies 



fh 



|A max -(A a + l)| < K-f^y) (38) 
where A a = max{7o,7i} and = min{7o,7i}. f(x) and 5 are defined as 

1 r ; n —\ll)T I 2 

f{x) == - [l-x-^{l-xf-^5\x\\ , 5 4 M L S '^ R ^ - « 1, xG (-oo,l-2 7 _]u[l+27+,+oo) 

with 7± = VS 2 +5±5. Furthermore, 0< f(x) <max{5, 7_} <C 1 when x G (— oo,l — 2^J\. 

Proof: The key idea here is to design the appropriate similarity transform for M in (31) and apply 
Lemma 2. Our detailed proof can be found in Appendix A. □ 

Now we are ready to derive the expression of G(SNR). We first consider the trivial case of — 1< 
7i<0. 2 Since 7d<---<7i» we can immediately have from (37) that O<V\s(A max )<^||-0T|| 2 ^l and 
0<^/(A max )< 7 ^||^ T || 2 <l- Therefore, (36) becomes G(SNR)«1 for all SNR. 

Next, we are going to discuss the case of 71 >0 in two separate cases. 

1) 7o > 7i •' Then A a = 7o, A;, = 71 and A max ~7o + l- So long as SNR is slightly larger than SNRt 
such that 7o>7i/(l— \/cr 2 / L\\iI>t\\)~7i, we have the following approximation for ips' 

2 D 2 D r^? 

u 7o ,2^ 7o <r v^i/ i2^./ cr 



^(A max ) - VE— i^i 2 < — < VtII^h « L 

L ^7o~7i 7o^i L f-' V L 



Moreover, we approximate tpj by the following bound, which is asymptotically tight with respect to ^a. 3 

\7o-7i/ V 70-71/ 



Substituting the above two approximations together with (34) and (20) into (36), we can finally have 

Pi + l 

Pil [1 - SNR T0 /SNRf + 1 
where Gy = I /(Pi + 1) and Pj is defined as 



G(SNR) = — OKir , - G v 



Pi± 



(N - ~ a 



(N - - 71 



l 2 ^2 D 

5^(7^ + 1)1^1 



L 

i=i 



2 Since 7; is the eigenvalue of (A;$AAf,Ri), 7»+l is for (cr^aoa^+Qs, Ri), a positive definite pair. Hence 7; + l>0. 
3 In fact, this approximation is precise enough when 70 is reasonably larger than 71, say, 70 > 271. 
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However, this expression is difficult to evaluate. Instead, we can compute Gjj first and have Pj=l/Gu—l- 
Noticing that Gjj is the output SINR as SNR=oo, we can ignore the term Aj^^Af on the left hand 
side of (23). Hence, w 00 =/iRj 1 ao= / u / Qj 1 ao, and by (17) and (14), we have its expression in (19). 

2) 70 < 71 : Then A a = 7i, Xb = Jo and A max ~7i+1- As a result, the term corresponding to 71 would 
dominate ijjg an d ipi in (37). To further evaluate these two terms, we need the bound in (38) to measure 



how close A max is to 71 + 1, namely, we use it to evaluate A D 



(71 + 1) in tp s and ipi: 



IpS (Amax)! 
■0/(Amax) ; 



71 



L A max -(7i+l) 
2 (7i+l)7i 2 



a 



1 



L 



IV'tJ 



(7 



L [A rj 



a 2 71 + I 



f SNR+l") 



71 



/ 70 



• 9 



L 



f 



71 + 1 



[^(Amax)] 2 



where g(x) = \ [l—x+\/ ( 1— x) 2 — 4<5x] . In fact, we can have the approximation that g(x) « 1 — x, because 

|(l-x)- 5 (x)| = |/(x)|«l 
Substituting the above expressions of ^5 and 7/7 as well as the approximation of g(x) into (3), we get 



G(SNR> 



N 

L/3SNR+7V 



1 2r 



71 
71-70 



^fSNR 



7i 



71+I 

xIVtJ 2 ^fSNR+l ' V7i-7o. 



N 

L/3SMR+N 



-1 2 



7i 
7i-7o 



1 2 



where G L = xlV"rJ 2 /(7i+l) is the output SINR when SNR=0 (-oodB). The last approximation holds 
when 7o<(l— \J q j^\i>T 1 D71, i-e. 70 is slightly smaller than 71, which implies 7 2 /(7i^yo) 2 <l/T~ 1 ^tJ 2 < 
(7i + IVxIV'tJ 2 - Furthermore, because ^ Tl | 2 <l, ^SNR/(^SNR+l)<K( 7 i + l)/^ Tl \ 2 . 
There is, however, no easier way to evaluate Gl than its definition. Fortunately, we are not interested in 
its specific values but the threshold SNRji that G reaches this value, as discussed in Sec. III-A. Noticing 
that SNRjq = 00 when ( N — (3)/f3 < 71 (c.f. (20)), we have following fact regarding 71 and SNRjo, 

N-f3 s 



= L(N - /3)SNR T0 . 
71 L/3SNR T0 + N ' 



Combining the above expression and (34), G(SNR) can be reduced to 

l + ^o 



G(SNR) 



- SNR/SNR T o + K (^SNR + l) 



G L 



(39) 



with Kq defined in (21). Finally, the following lemma provides an easier way to compute SNRjo from 
71. The proof can also be found in Appendix A. 

Lemma 5: 71, ... ,jd are approximate to all the nonzero generalized eigenvalues of (Qs—Qx, Qx) 
(padded up to D with zeros if not enough), and is almost independent of SNR. 
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Input SNR(dB) SNR (dB) 



(a) The curves of 70+I, 71+I and A max vs. SNR. Q3 = 0) (b) 70 + I against SNR with different /3. 

Fig. 3. The curve of 70 + 1 and 71 + 1 against SNR. The intersection of them is the empirical threshold SNRtq. 



IV. Discussion of the Threshold 

From Theorem 1 in Sec. III-A, we know that the empirical threshold SNRjo is a key parameter for 
MPB. By the discussion in Sec. III-C, it is the intersection of the curves 70+I and 71+I, as shown in 
Fig. 3(a). Therefore, we need to investigate these two parameters to gain deeper insight of SNR-r - 

A. General Results of 70 and 71 

70 has a simple expression of (34), from which we know that the parameter f3 is critical. Fig. 3(b) 
shows the curve of 70 + 1 with different (3. We see that, if /3 = 0, i.e. Hj and are orthogonal, then 
7o=LSNR is unbounded as SNR goes to infinity. Otherwise, there would be a limiting value of (N—j3)/0 
so that 70 might never exceed 71 and SNRjo=+oo. Therefore, 0=0 is the best choice for 70. 

71 is another critical parameter that determines SNRjo- Fig. 3 shows that, as 71 increases, its intersec- 
tion with 70 moves rightward and SNR-rn increases. Therefore, knowing how to control 71 is important 
in designing MPB. By Lemma 5, 71 is the largest nonzero generalized eigenvalue of (Q5— Qx, Qx)- (It 
is zero if there are less than D nonzero eigenvalues and all of them are negative.) To solve it directly, 
we need the roots of a polynomial eigen-equation of order L, which has no general formula when L>4. 
Instead, we resort to matrix perturbation theory again to derive a bound for it. We first notice that the 
eigen-decomposition of (Q5 — Qx, Qx) is equivalent to (Qs,Qx) and their eigenvalues only differ by 
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one. By (6) and (7), the eigenvalue of (Qs, Qx) is further equivalent to that of (Ys+I, Yx+I), where 



Y s = ^A 7 * 5 Af = INR • A/* 5o A? , Y x = ^A 7 * x Af = INR • A 7 *x Af . 



(40) 



We term (Ys, Yr) as noise free pair, since it can be viewed as the covariance matrix pair of MPB 
without noise. Our strategy here is to view (Ys + I, Yr + I) as a perturbed version of (Ys, Yr) and 
apply the results in matrix perturbation theory to derive a bound for 71+I. Before we proceed on, we 
cite a more general definition of the generalized eigenvalue of a matrix pair [32]. 

Definition 2 (Generalized eigenvalue): The generalized eigenvalue of a matrix pair (A, B) is a one 
dimensional space, denoted as (v,p) = {r ■ [v p] : r € C}, where [v p] is a 1 x 2 vector satisfying 

p • Ax = v • Bx, 

with x being its eigenvector. If B is nonsingular, then \=vjp becomes the conventional definition. 

Comparing to the conventional definition of generalized eigenvalues, this one includes the special case 
of B being singular so that ^=0 and v^O, namely, A=+oo. We will see its importance later. Besides, 
we also need the following definition and lemma from [32, pp.3 15-3 16]. 

Definition 3 (Definite pair): If a matrix pair (A, B) consists of two Hermitian matrices and 



C(A,B) = min i/(x H Ax)' + (x^Bxf > 0, 

||x||=l V 

then (A,B) is called definite pair, and C(A, B) is its Crawford number. 

Theorem 2 (Weyl-Lidskii type): Assume (A, B) is a definite pair and (A+E,B+F) is its perturbed 
version. Let (aj,/3j) and (i=l, . . . ,ri), be their ordered generalized eigenvalues, respectively. If 



+ ||F||| < C(A,B), (41) 
where || • \\s denotes spectral norm of a matrix, then (A, B) is a definite pair, and 



\ai02-Pia 2 \ < V I|E|I ^ + I|F| 



VI«il 2 + l/3ilVl5 2 | 2 + |/3 2 | 2 C(A ' B) 

We want to apply Theorem 2 to (Ys+I, Yx+I) and (Ys, Yx) to derive a bound for its eigenvalues. 
However, (Ys,Yx) is not a definite pair, for the null spaces of Ys and Yx may have nontrivial 
interset (larger than {0}). Therefore, we need to transform (Ys, Yx) into a definite pair. Let Ms and 
Mi be the null spaces of Ys and Yx, respectively, and Mo = Ms fl-^r- Let Eo be a matrix whose 
columns are the orthonormal basis of M^, with E^Eo = I. Then, by the determinant identity in (29), 
(Ys+I, Yx+I) has the same generalized eigenvalues as (E^YsEq+I, E^YxEq + I), except for the 
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multiplicity of ones. Similarly, (E^Y^Eo, E^YjEo) has the same eigenvalue as (Ys, Yz) does except 
for the multiplicity of zero. Therefore, instead, we can apply Theorem 2 to (E^Y^Eo, E^YjEo) and 
(E^Y 5 E +I,E^Y 5 E +I) to derive the bound. The Crawford number of (E^Y^Eq, E^Y X E ) is 



Cy= min J(w H Y s w) z +(w H Y x w) > 0, (42) 
weA/i 1 , || w||=i v 

which means (Eq^Y^Eo, E^YjEo) is a positive pair. Before applying Theorem 2, we first analyze the 
dependency of Cy on the interference power. Substituting (40) into (42), we have Cy = INR • Cy , where 

C Ya = min A/(w^A/$ 5o Afw) 2 +(w H A/*x Afw) 2 . 
weA^, ||w|j=i 

Since &s an d $>i are independent of INR (c.f. Sec. ITB), Cy a is also independent of INR. Therefore, Cy 
is proportional to INR. Now, we are ready to use Theorem 2, and we only consider the case of INR being 
large. Let (uq, /io) and (v, fi) be the corresponding generalized eigenvalues of (E^Y^Eo, E^YjEo) and 
(E^Y 5 E +I,E^YxE +I), respectively. Then, for large INR, (41) is satisfied and 

\ufi -fiu \ < ^2l.J_ (43) 



VlH 2 +HM^TO^F °y a INR ' 

We now derive the bound for 71 + 1 in two separate cases. 

1) Mi ^ Ms-' There is a nonzero wo such that Y5W0/O and Yjwo = 0. By Definition 2, (fo,0) 
(Wo^O) is a generalized eigenvalue of (Y 5 , Y x ) and (E^YsE , E^YjE ). This means the noise free 
pair (Ys, Yz) has an infinite generalized eigenvalue. As a result, (43) becomes 



VM 2 +M 2 C Yo INR n V 2 V2 

This means (E^Y^Eo + 1, Eq^Y^Eo + 1) would always have an eigenvalue that satisfies the above 
inequality. Since (Qs, Qz) has the same eigenvalue except for the multiplicity of ones, 71 + 1 satisfies 

7l + l>^|lNR, (44) 

which gives a lower bound for 71 +1. We can see that it goes up unboundedly as INR increases. Hence, 
by (20) and (22), the threshold SNR also increases unboundedly with INR. 

2) Mi C Ms-' Then, Mo =Mi and E^ YjE is nonsingular. As a result, all eigenvalues of (E^ Y5E0, E^YjEo) 
and (Ys,Yx) are bounded. When INR is large, (43) becomes 

|A " A0 ' <%-i4 =► |A-A |<%(1 + A 2 ). 1 



^/TTaVT+AI V2 INR ' Ul V2 y u/ INR 

where Ao = ^o/Vo is the largest eigenvalue of (Eq^Y^Eo, Eq^Y^Eo) and we used A« Ao in the second 
inequality. Since (Qs, Qz) has the same eigenvalue as (E^Y^Eq + 1, E^Y^Eq + 1) except for the 
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multiplicity of ones, 71 + 1 is bounded around the largest eigenvalue of (Ys, Yx) or one. Furthermore, 
from (40), we know that Ao is independent of I NR. Therefore, as all eigenvalues of noise free pair is 
finite, 71 + I is bounded and independent of I NR. 

B. Two Typical Scenarios of 71 

From the previous part, we know that whether the threshold of MPB is unbounded is determined 
by the existence of infinite eigenvalue of the noise free pair. Now, we discuss two typical classes of 
interferences. In the first case, all eigenvalues are finite, and in the second case, there may be infinite 
eigenvalues. 

1) Directional White Noise: By uncorrected directional white noise, we mean an interference that 
arrives from a specific direction, with its samples in time domain being uncorrected. This means the 
entries of S/(fc) are independent identically distributed random variables with zero mean and unit variance. 
(In section II-A, we have already normalized the interference power in Sj(fc).) Then, by (8), we can 
immediately have $>s a = &z , and all the eigenvalue of the noise free pair (Ys, Yx) are ones. Thus, in 
this case, 71 + 1 is bounded and is independent of I NR. In fact, we can further have Qs = Qx according 
to (6) and (7). Therefore, 71 = 0, SNR TO = and G(SNR) = 1. 

2) Directional Periodical Interference: If an interference has periodical structure in time domain and 
arrives from certain direction, then we term it as directional periodical interference. By periodical, we 
mean the interfer is periodic with respect to the projection basis, i.e. 

Sj(fc) = e^HfSj, hfSj(fc) = e^hf Sj. (45) 

A stronger condition would be S/(fc) = Sj. However, (45) is good enough for our analysis. Now, we will 
discuss the existence of infinite generalized eigenvalue in the noise free pair (Ys, Yx). To do this, we 
need to consider the relationship between their null spaces. 

In practice, Rs and Rx are estimated from sample average. Therefore, replacing the expectation in 
(8) by sample average and using (45), we can get 

*s = <sj h^s^f, *x = In^sfHjH^sjny 2 . 

Let n = S/fiy 2 Aj. Then, by (40), Ys and Yx in this case can be expressed as 

Ys = INR ■ (n^hshfn)* , Yx = INR- (n^HxHfn)*. (46) 

Therefore, the noise free pair (Ys, Yx) has the same eig envalue as (Il^hshf II, ^-II^HxH^II) does, 
and we only have to check the null spaces of later matrix pair. Let Vi = 1Z(I1) be the range space of II. 
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By definition, it is also the range space of Sj, namely the space spanned by all interference waveforms 
in one period. Then, we can express the null spaces of n^h^h^n and II^HiH^n as 

M^=M (u H h s h^u) = {nt x: xe5 1 nv J }®JV(n) 
N' z = TV (n H HiHj n) = {n+x: x e^n v z \ eAA(n), 

where S = TZ(hs), Z=lZ(Hz) as defined in Sec. II-B, and f denotes Moore-Penrose pseudoinverse [34]. 
To see the relationship between M' s and Wj, we only need to check the relationship between S 1 - n Vi 
and X L n Vi which can be expressed as 



5 1 nV/ = {x: Bf x = 0}, I 1 nV i = {x: B^x = 0}, 



(47) 



where 



B x 



H 



H 7 



B, 



with H v _l being the matrix whose columns are the orthonormal basis of Vj~ . Then, by (47), whether the 
eigenvalues of (Y^, Yj) are finite is equivalent to the validity of TZ{Bx}^JZ{'Bs}- Define 



I -HylHj 

I 



I -H^h 5 
I 



Right-multiplying Bj and B^ by Tb z and Tb s , respectively, we can have 



Bi • Tb z 



B5 • TB ^ 



where Py 7 ^H-^Hy is the projection matrix onto V/. Since Tb x and Tb s are nonsingular, the above 
two matrices have the same range spaces as B^ and Bj do, namely we can have 

TZ{B X } = Vf TZ{P Vl U x }, TZ{B S } = Vf H{V Vl \v s } 

Therefore, whether the generalized eigenvalues of the noise free pair are bounded is finally equivalent to 
whether 



ft{p Vj H x }Dft{Pv c h lS }. 



(48) 



Fig. 4 shows the geometrical interpretation of (48). To see if 71 is bounded, we can project h,s and 
all the columns of Hj onto V/. If the former projection lies in the space spanned by the later ones, 
then 71 is bounded. Otherwise, there is an infinite eigenvalue in the noise free pair and 71 will increase 
unboundedly with I NR. Since Hj of the Maximin scheme and PAPC scheme has only one column, this 
means the projections of hs and Hj should be on the same line, which is hardly valid. Thus, these two 
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One Interference Channel Two Interference 
Co (Maximin / PAPC) c Channels 




V, = Tl{s,} is the subspace 
spanned by the waveforms of the interferers in one period 



Fig. 4. Geometrical interpretation of (48). and h^ 2 ' denote the columns of Hi. For one interference channel, Hi = hi. 

approaches are sensitive to directional periodical interferences, as we will see in later simulation results. 
Furthermore, Fig. 4 also shows the case of multiple interference channels, i.e. the dimension of X is 
larger than one, then 71 is more likely to be bounded. This issue is out of the scope of this paper and 
will be discussed in in [36]. 

V. Simulation Results 

In this section, we simulate various scenarios and compare them to our theoretical results. The matrix 
pair beamformers implemented include the Maximin algorithm and PAPC algorithm, which have many 
kinds of adaptive algorithms. However, we are only interested in their steady state performance. The 
reason for this is that, if the methods suffer from poor steady state performance, then it is meaningless to 
investigate their adaptive algorithms. Therefore, we directly calculate their weight vectors by performing 
generalized eigen-decompositions on the estimated matrix pairs. The interference signals encountered in 
the simulation include BPSK signal, tones, periodical noise and multiuser interference. The first one is 
the directional white noise and the others are directional periodical interferences. 

In all cases, we consider a uniform linear array (ULA) with eight isotropic antennas (L = 8) spaced 
half a wavelength apart. For each user, a 100 kbps DPSK signal is randomly generated and spreaded by 
a distinct 31 -chip Gold sequence in each simulation trial. Then it is modulated onto a 1 GHz carrier to 
form a RF signal with bandwidth 3.1 MHz. In all the simulations, we assume the SOI arrives from 0° 
and the interferers have equal power. 

Fig. 5(a)-Fig. 5(d) show the operating curve G(SNR) of the Maximin and PAPC beamformers in the 
four scenarios. The simulated values are obtained by using the simulated data and the theoretical ones 
are computed by the piecewise function G(SNR) in (18). To eliminate the randomness caused by finite 
sample effects, K = 10 6 symbols are simulated for each SNR and INR in every experiment. 
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(a) G(SNR) of Maximin and PAPC under BPSK jammers, (b) G(SNR) of Maximin and PAPC under periodical noise. 
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(c) G(SNR) of Maximin and PAPC under tone jammers, (d) G(SNR) of Maximin and PAPC under multipath MAI. 
Fig. 5. Simulated and theoretical operating curve G(SNR) under different kinds of interference for Maximin and PAPC. 



In Fig. 5(a), the interferers are three uncorrelated BPSK modulated signals, i.e. random sequences of 
±1. The rates of the interferers are all 3.1Mbps, which is the same as the chip rate of the SOI and 
covers its whole bandwidth. Therefore, they belong to the type of directional white noise. We can see 
that G(SNR) w 1 and there is no threshold effect, which is consistent with the analysis in Section IV-B1. 

Fig. 5(b)-Fig. 5(d) show the results for three types of directional periodical interferences: periodical 
noise, tones and multiple access interference (MAI) with multipath. In Fig. 5(b), we consider two 
periodical noises, which arrives from 30° and —40°, respectively, and each of them is generated by 
repeating a segment of Gaussian white noise over times, with the repeating frequency being 100 kHz 
(same as that of SOI's symbol rate). The frequency offsets of the tones with respect to the carrier in Fig. 
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5(c) are 100kHz, -300kHz, OHz, 400kHz and -100kHz, respectively, and their DOAs are 30°, -50°, 
—20°, 19° and 45°. In Fig. 5(d), there is one incident MAI signal with three -ray multipath delays of 3 
chips, 5 chips and 4 chips arriving from 30°, —20° and —50°, respectively. We can see that the theoretical 
values of G(SNR) match with the simulated ones very well. 

To verify the validity of the approximation given by Lemma 4, in Fig. 6(a) and Fig. 6(b), we also 
show the curves of 70 + 1, 71 + 1, and A max . We can see that max{7o + 1,71 + 1} can be an excellent 
approximation for A max . In Fig. 7(a) and Fig. 7(b), we show the array patterns that are below and above 
the threshold. The results in these four figures are simulated under two periodical noises, with the same 
parameter as Fig. 5(b). The case for the tone jammers and the multiple access interference cases are quite 
similar, and is thus omitted. 

For the Maximin algorithm, all the curves of G(SNR) have failure area, threshold area and operating 
area, which are consistent with the typical curve in Fig. 2(a). We have also marked the predicted thresholds 
of SNRji and SNRj2, computed by (22), in the figures as well, which confirm our theoretical calculations. 
Furthermore, they also show that the threshold SNR would increase with the interference INR. This is 
because 71 + 1 (c.f. Fig. 6(a)) would moving upward as INR increases, making the intersecting point of 
70 + 1 and 71 + I move rightwards. This is consistent with our claim in Section IV-B2 that 71 would 
increase unboundedly with the INR. Fig. 7(a) and Fig. 7(b) show its array pattern below and above 
threshold, respectively. We can see that the mainlobe of the beamformer would mistakenly point to the 
interferers once the SNR is below the threshold, as we have predicted in Section III-C. 

For the PAPC algorithm, Fig. 5(b)-Fig. 5(d) show that there are only failure areas. This is because, 
no matter how large SNR is, 70+I will never exceeds 71+I, as shown by Fig. 6(b). Therefore, the 
threshold of PAPC in this case is infinity. The array pattern shown in Fig. 7(a) also indicates that its 
mainlobe has pointed to the interferers. We have stated in section III-C that the presence of SOI in Rj 
will make the beamformer mistakenly null the SOI. The curves of PAPC in Fig. 5(b)-Fig. 5(d) confirms 
this observation, and Fig. 7(b) shows that PAPC beamformer has a deep null in the direction of the SOI. 

In summary, the conventional MPB like Maximin and PAPC work well in the presence of directional 
white noise, even when the INR is large. However, they are very vulnerable to multiple directional 
periodical interferers with repeating structures in the time domain and arrive from different directions. 

VI. Conclusions 

Matrix pair beamformer (MPB) is a general framework we proposed to model a class of blind 
beamformers that exploit the temporal signature of the signal of interest (SOI). It has the advantages 
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(a) Maximin scheme 
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Fig. 6. The largest and second largest generalized eigenvalues of the Maximin scheme vs SNR under two periodical noises. 
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(a) Case I: below threshold (SNR=-10.1dB, INR=30dB) (b) Case II: above threshold (SNR=40.9dB, INR=30dB) 



Fig. 7. The array patterns corresponding to the Maximin and PAPC under two periodical noises. 



only relying on the second order statistics to achieve blind processing. In this paper, we have analyzed 
the mechanism of MPB with matrix mismatch, and showed how it worked "blindly". We have discovered 
that there is a threshold effect in MPB, i.e. the beamformer would fail completely if the SOI's input SNR 
is below that threshold. Meanwhile, its normalized output SINR has been derived as the performance 
measure, and the threshold SNR has also been predicted. We have also observed that the existence of 
infinite generalized eigenvalue in what is called noise free pair makes the threshold increase unboundedly 
with the interference power. This is highly probable when there are multiple periodical interferers. All 
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our theoretical analysis matches with the simulation results very well. 

Our analysis indicates that the conventional MPB is very vulnerable to multiple periodical interferers. 
Moreover, it also implies the importance of choosing the appropriate projection space for the interference 
channel. And we will address this issue in another paper [36]. 

Appendix A 
Approximation of A max and 7, 

A. Approximation of A max 

Since similarity transform does not change the eigenvalues, we can apply it to M before using Lemma 
2. As we will see later, the bounds derived in this way can be surprisingly tight. Though any transform 
matrix can be used, we prefer the following diagonal matrix that can preserve the diagonal terms of M 

F = diag{l,/ 1 ,/ 2 ,...,/ I3 }, (ft >0, i = l,2,...,D) 

Applying it to M in (31) and using Lemma 2, we can have the following Gerschgorin disks for FMF" 1 : 

G l = {\: |A-( 7i + l)| <Ri(F)}, i = 0,1,2,..., D, 

where the radii of the disks are 

D 1 

R^(F) = l0 J^-J2\^ i \j (49) 
i=i 

, . VxIV'tJ ■sgn(7 i ) 

Ri(F) = "fi^-^a fi, l<i<D, (50) 

; u ^SNR + 1 3 ~ ~ 

with tpT z being the ith component of Vt- Now, we are going to optimize /i,...,/d- To derive an 

effective bound for A max , we should ensure the rightmost disk is separated from the others (c.f. Lemma 

2). Therefore, our criterion for finding the optimal F is minimizing the radius of the rightmost Gerschgorin 

disk subject to the constraint that it is separated from all the remaining ones. Depending on 70 + 1 > 71 + 1 

or 70 + 1 < 71 + 1, the rightmost disk might be Go or G\, and we now discuss them separately. 

If 70 + 1 > 71 + 1, then Go is the rightmost disk, and we can formulate the optimization of F as 

r~2~ D 1 

min R (F)= l0 J^-J2\^\j (51) 

i=i 

s.t. ji + l + Ri(F) < 7o + l- J Ro(F), l<i<D. (52) 

To minimize Rq(F), we want f±, . . . , fjj to be as large as possible. But this will increase the radius of 
Gi, . . . ,Gd, and to avoid connecting Go, they cannot be too large. In fact, /1, . . . , /d have different 
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importance in this tradeoff. Since 70 > ... > 713, G2, ■ ■ ■ ,Gd are farther away from Go than G±. 
Therefore, /2, • ■ ■ , /d can be reasonably larger than fi when keeping separate from Gq. Thus, to simplify 
analysis, we can ignore all the terms except (V'tJ/Zi in -Ro(F), i.e. Ro(F) « 7o^/^"^~ m (51) and 
(52), and only the constraint for i = 1 is effective in (52). Then, the problem can be reduced to 

[o* 1 
mm Ro(F) « 7 y — \ip Tl \ ~r 



v xIV'tJ -sgn(7i) /^2~ 
S>t - 71 ^SNR + 1 ^ " (70 "7i)/i + 70^1^1 < 0. 

This is a simple convex optimization problem, feasible when ^ G (—00, 1 + 2(5— -^/(l + 2(5) 2 — 1], 
where 5 = ^|^Ti| 2 /(t^SNR + 1) < 1. Thus, we can easily solve the above optimization and get 

|A max -( 70 + l)| < 7 o-/(^) 

where f(x) = \[\-x- y/ (1- x) 2 - 48\x\] with x G (—00, 1 - + 25} U [1 + 2 Vd+5 2 ~+ 25, +00) . 

Furthermore, by taking the derivative of f(x), we can easily get 0</(x)<max{5, \/5 2 +5—5}<^il when 
x € (-oo,l-2V*+5 J +2(5]. 

If 70 + 1 < 71 + 1, then Gi is the rightmost disk and 71 > 0. The optimization problem is 



mto = 7i MsnrTT <53) 

s.t. 70 + I + Ro(F) < 7l + l -i?i(F), (54) 
7i + l + i2i(F)<7i + l-i2i(F), 2<i<D. (55) 

According to (49) and (50), /1 should be as small as possible to minimize -Ri(F) in (53), (54) and 
(55). However, this will increase Rq{F) in (54), making Go connect with G\. Therefore, f\ cannot be 
arbitrarily small. On the other hand, fz, . . ■ , /d should be as large as possible to reduce Rq{F) in (54) 
while keeping (55) valid. As a result, \ipT y \l h is still the dominant term in i?o(F) (c.f. (49)) and the key 
point here remains the tradeoff between Ri(F) and Ro(F). In other words, the optimization becomes 

min i?i(F)=7i T -5 U 

u ' ri M S NR + l 



s -'- "1 /. ,g NR ,/i -(71 -7())./'i +7()\ , ' < »■ 
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The feasible region for this convex optimization is ^ G [1 + 25+ \/ (1 + 25) 2 — 1, +oo), with 5 defined 
in the previous case. By solving it, we can finally get the bound for A max as 

|A max - (71 + 1)| < 71 ■ / 
where f(x) is the same as in the previous case. 




B. Approximation of 7$ 

Then, we discuss the approximation of 7^. By substituting (5) into (28) and using matrix inversion 
lemma, we can have the equivalent equation of (28) as 

erf Af Q-^o • a^Q^Aj 



det(Af R^A/)" 1 • det <^ AI 



AfQ^Aj 



*A > = (56) 



l + afa^Q^ao 

To further reduce the above expression, we need factorizations of AfQ^Aj and <1>a like (26). Let 



o\ = in (26). Then, Rj = Qj, and (26) become 



'T 
Ml 



T^* A T = r , T?(Af Q^A,)" 1 ^ = I. 



where To and To are the counterparts of T and Y, respectively. Substituting them into (56), we have 



det(AfR~ 1 A / )- 1 - det <^ AI 



1 + of a^Q^ao 







where ipT = T ~ 1 A^Qj 1 ao is the counterpart of ip? when a\ a = 0. By letting SNR = in Lemma 
1, Rj becomes Qj and we can have a^Q^ao ~ L/a 2 and H'j/'ToII 2 = Lko/ct 2 <C L/cr 2 . Therefore, 
by the similar argument of A max using Gerschgorin theorem, we can have 7, « A^o, where A^o is 
the ith diagonal term of To, namely, the ith eigenvalue of (3>a, (A^Qj 1 A/) -1 ). In fact, there is a 
correspondence between A^o and the generalized eigenvalue of (Q5 — Qj, Qj). By (29), 

det {A(Af Q^A/)" 1 - * A } = \ D ~ L ■ det(Af Q^A/)" 1 det Q^ 1 • det{AQ x - A/$ A Af }. 

Since Q5 — Qj = A;$AAf, the above expression implies that (<J>a, (A^Qj 1 A/) -1 ) has the same 
eigenvalues as (Q5 — Qj, Qj) except for multiplicity of zeros. Thus, we can estimate 7, like this: 1) 
take out all nonzero eigenvalues of (Q5 — Qj, Qx), 2) pad them up to D eigenvalues with zeros, 3) 
order them decreasingly to get Ai ; o, • • • , Ad,o> an d 4) let ji « A^q. 
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Appendix B 
Proof of Lemma 3 

Substituting the definitions of ys(k) and yi(k) (c.f. (10)), (6) and (14) into (17) , we can have 

fir c Nm = |w g a | 2 1 __ „ |w g a | 2 /A 

{ ' w^Aj^Afw + ^Hwfa^Q^ao ~ w^A/^Af w + a 2 ||w|| 2 / a 2 ' yJ,} 

where we used a^Q iS 1 ao~L/<r 2 , derived by replacing Rj with and letting o 2 s =0 in Lemma 1. 

To further derive G, we first need the expression for ||w|| 2 . The key trick is to recognize that w G 

71(A) so that its projection onto H(A) equals itself, where IZ(-) is the range space of a matrix and 

A = [ao A/]. This can be proved by substituting (7) into (5), applying matrix inversion lemma and 

pluging it into (25). Furthermore, IZ(Ai) is a subspace of 11(A) with the dimension lower by one. 

Thus, 1Z(A)=1Z(A)(BlZ(bo), where © denotes the direct sum and bo is the unit vector in 1Z(A) that is 

orthogonal to 1Z(Aj). Then, the projection matrix of 1Z(A) can be written as 

P A = A 7 (Af A,)- 1 Af + b b^. 

This together with Paw = w and P 2 A = Pa yields 

||w|| 2 = w^Paw = w H A/(Af A 7 ) _1 Af w + |w H b | 2 (58) 

Assumption 1 in Sec. III-A implies that ao is almost orthogonal to 1Z(Aj). Thus, it is nearly aligned 
with bo, and intuitively, jw^bol ~ |w^ao|, where ao = ao/||ao|| is the unit vector of ao- To prove it, 
let IZq = spanjao, bo} and {h>o,bi} be the orthonormal basis of IZq. Its projection matrix becomes 
Ptto=bobjf+bibf , which satisfies P7e o a =a and P7e o b =b . Define wo=Ptc w/||P^ w||. Then, 



|w H b | 2 = |w^a | 2 + 



w H bo| 2 — Iw^aol 



\w H a \ 2 + 



| P^ b 1 2 -| w H Pn„ Pn a 1 2 



= |w^a | 2 +||P7e o w|! 



|w^b | 2 -|w^b • b^ao+w^bi • bf a | 2 

= |w^a | 2 +||P TCo w|| 2 • cos fa • cos(0i + 20/) (59) 

where 4>j is the angle between ao and 1Z(Aj), and cf>\ is the angle between ao and bo- And the derivation 
of the last step in (59) involves some simple trigonometry identities like product-to-sum formula. Let 
tf 7 =Af A//L. Combining (58), (59), || Ptc w||<||w|| and | cos^>/|<Cl (c.f. Assumption 1), we can have 



Iw|l 2 = 



f la^wl 2 + A/tf 7 1 Af w 1 n „ l9 „ a t , a „ n 

l-^^-cos</> 7 -co S (</> 1 + 20/) L 



| w 

2 



where we also used ||ao|| —L in Assumption 2. Substitute the above expression back to (57), and we get 

G(SNR)« — rj |w ^ o|2 „ , - - . (60) 

v ' w^A / [^$ lS +^7 1 ]Afw+|w^a | 2 
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Next, we are going to evaluate w^ao and w^A^^^Ftj 1 ]A^w. By the expression of w in (25), 

D D 



a^w = a^R^ao + J] 



i=i 



At: 



-|^| 2 = a ^R x 1 a -||^ T || 2 +^ 



A T1 



-1 



L (LP 



^SNR + 1 



a 2 \N 



c -(7i + l) 

1 + (^SNR I ) r,(A ma ,-) 



i=i 



A, 



"(7i + l 



(61) 



where we used Lemma 1 in the approximation and ^s(A max ) is defined as (37) in Lemma 3. Before the 
derivation of w H A/[^<l> l s+\l/7 1 ] A^w, we first cite the following identity from Lemma 2 in [35]. 



L/3SNR • ^t*Pt 
I_ L/3(l-p )SNR+iV 



T" 1 . 



Then, combining the above expression together with &a = &s— ^x, (26) and A e =A/(T 1 ) H we can 
have 



L 



L 



Afw = -^w H A I (T~ 1 ) H 



NL/a 2 



L/3SNR • ^ t iI>t 
LP{l- Pa )SUR+N 

■^/(Amax) 



i+r- 



T _1 Afw 



L/3SNR 2 

V's(Amax) 



L/3SNR + iV 



(62) 



L/3SNR + iV 

where ^/(Amax) is also given by (37) in Lemma 3, and po<^l is a small number independent of SNR 
Finally, substituting (61) and (62) into (60), we can derive G(SNR) in (36) immediately. 
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Abstract 

We give some key results for the paper "Performance Analysis of the Matrix Pair Beamformer with 
Matrix Mismatch". 

Index Terms 

Adaptive beamforming, Matrix pair beamformer, MPB, Generalized eigen-decomposition, Matrix 
mismatch 

I. Introduction 

We give some key results for the paper "Performance Analysis of the Matrix Pair Beamformer with 
Matrix Mismatch". 

II. Problem Formulation 

The analysis of the performance of the matrix pair beamformer requires the expression of A // Rj 1 A, 
where 

R x 4 a| o a a^ + A 7 * x Af + a 2 I, (1) 



a ai • • • a D 



a A/ 
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Furthermore, we base our discussion on the following assumptions. 

Assumption 1: The spacing of DOA between any two signals is large enough (greater than a mainlobe), 
so that {aj}^ is linearly independent and the projection of ao onto span{aj}^ =1 is much less than ||ao||. 

Assumption 2: The steering vector of all signals are normalized so that ||aj|| 2 = L, (i = 0, 1, . . . , D). 



III. Some Fundamental Results 

The main results are summarized in the following two Lemmas. 
Lemma 1: The expression of A H YL I 1 A is given by 

a^R^ao a^R^A/ 



A^R^A 



where 



a^R 2 - 1 ao - 



AfR-W 



AfR^Ar 



L 



AfR^ao A^Rj 1 A/ 



° 2 t#(1-0SNR + 1 



L 



1 



a2 ¥(1-0SNR + 1 



(7 



Mqmr 
f (1-0SNR+1 J 7 



(2) 
(3) 

(4) 



with ip 4 ^r 7 A A^i. ^ = Naj .j aSo . B and ^ ^ defined as 



£ - Po - «0- 



po — ip H ty x 1 il> is the projection power of the unit vector of ao onto spanjao, 
smaller than 1; kq = ip H ^J 1 E i ^J 1 ip. Furthermore, £ and kq satisfy 

< «o < p < 1, < £ < po < 1 

Proof: According to (1), Rx can be written as 

R x = ASA^ + a 2 I 

where A and S are defined as 

a A/ 



, &£)} and is thus much 



H 



a* 

-La 







(5) 



(6) 



(7) 
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According to matrix inversion lemma [1, pp.124 (3.8.3)], we can derive the R^ 1 as 

" L 



R-7 



L jl 



L l -1? K V' S 



*S + I 



(8) 



where \£ = A H A/L. Let tp = Afa /L and \£/ = Af Aj/L, then * can be written in the following 
block form 

"l t/>" 



(9) 



By (8), we can write A H ~R X 1 A as 



A^R^A = ^ S * " 



(7 



(7^ 



(7 



#S + I 



(10) 



In order to further derive the expression for A H R X 1 A, we need to write [\P • LS/cr 2 + 1] 1 in the 
following block form 



L 



*S + I 



X\\ X12 
X21 X22 



(11) 



where x\\, X12, X21 and X22 are 1 x 1, 1 x D, D x 1 and D x D matrix, respectively. In order to derive 
their expression, we need the following equality [1, ppl23] 



(12) 



where S = B — RA X C is the Schur complement. Substituting (7) and (9) into (11) and using (12) as 
well as matrix inversion lemma, we can derive 



A 


C 


-1 


"a^+a^cs^ra- 1 -a^cs- 1 " 


R 


B 




-S-iRA- 1 S- 1 



arn = (1 - a) 



a 



xi2 = -(l-a)-^-^*zP 



x 2 i = —ol\ ■ Ftp 

X 22 = P + 0,4 • P^ H *XP 



(13) 
(14) 
(15) 
(16) 
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where 



L 



(7 



+ I 



a 



La 1 



a' 



+ 1 



Ql 



(7 



La 2 

-M 

a 2 



-(i-O + i 



(17) 
(18) 

(19) 
(20) 



It is easy to prove that a, a±, £ and P satisfy the following equalities 

ai = a(l + ai£) 



P H = I — ^$jPf 



which have already been used during the derivation of (13)— (16) and will be used to simplify the 
expression of A H H X 1 A. Substituting (9), (11), (13)— ( 16) into (10), after some lengthy but straightforward 



derivations, we can have 



x a 1 



(l-a)^P^ Ptf 7 -aiP^ H P H 



(21) 



Using (6), (17)-(20), SNR = a 2 s Ja 2 and /3 = Na^-Ja^ we can finally derive (2)-(4). 

The remaining problem is to prove £ = po — kq and k G [0,po], which can also naturally yield 
£ G [0, po]- According to (20) and the definition of H, £ can be written as 

= ^ H s _1 s*7V - v^tf 7^*7 V 



Po 



K 



Since 3>j is positive semi-definite, H 7 1 — ^ 1 = is also positive semi-definite, i.e. \t 7 1 < H 7 1 . 

Then */ > Hj, thus 



£ = V"*7 



*7V > o, 



which means po > K o > 0. By the definitions of ?/> and we have po 



1 n H 



a^A/(Af A/J-^jao, 



i.e. the projection power of ao / 1 1 a o 1 1 onto span{ai,a2, . . . ,sld}, which we have assumed to be much 
smaller than 1 in section II Thus po <C 1 and this completes the proof. □ 
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Lemma 2: Let i})t — 7^ (1 — Po)SNR + 1 ip T with ip T = A e R :r 1 ao, then 



L 



L_ g(i-po)SNR + l 
7? (1-OSNR + l 



-H 



a'- 



4?SNR • y~ , iH 



f (1- P0 )SNR+1 



(22) 
(23) 



where po, £ and kq are defined in Lemma 1, and are independent of SNR. 

Proof: Since A e =AjT~ H and Af R^ 1 A/=TT H (c.f. (26) in [2]), (4) is equivalent to 



M (1-0SNR+1 



A/ 



TT 



H 



(24) 



Right-multiplying it by ^ and using po = £ + K o> we can derive 



V = 



a 



2 M 

N 



(1-OSNR + l „ , 



£ f(l- P0 )SNR + l 



(25) 



Since A e = A/T H , then by (3) and above expression, we have 

L 1 



$ T = Af R-^o 



a2 t#(1-0SNR + 1 



T _1 S¥7 V = 

^(l-p )SNR + l 



Let ^ 



(l-p )SNR+l Vt» then i/'t = T // ^ / 1 i/»- Using (24), p =Ko+£ and the definition of k 



(Lemma 1), we have 



\\^ T f = t/> // fl' / 'tt"*, 'v = 



L ^(l-p )SNR + l 



cr 



_N_ 

2 L/3 



(1 -0SNR + 1 



K . 



According to Lemma 1, 0<£<Cl, 0<po<^3- and 0<t,/ko<^/po. Therefore, HV'rlNy ^ K o^y We then 
prove (23). Substituting (25) into (24) and using V>t=T h \E , 7V' we can have 

^1-0SNR+1 



H = T< 



L 



■1+ 



Mi 

N 



^SNR 

2 — V'tV't 



M ( i_ po) SNR+l 

Applying matrix inversion lemma [1] to the above equation and using po = K o+£ (Lemma 1), we can 



finally derive (23). 



□ 
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